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Several results are presented concerning the equation F(x) + Ax = B 
(with F(-) a "diagonal" nonlinear mapping of real Euclidean n-space E" 
into itself, and A a real n X n matrix) which plays a central role in the 
dc analysis of transistor networks. In particular, ice give necessary and 
sufficient conditions on A such that the equation possesses a unique solution 
x for each real n-vector B and each strictly monotone increasing F(-) that 
maps E n onto itself. 

There are several direct circuit-theoretic implications of the results. For 
example, we show that if the short-circuit admittance matrix G of the linear 
portion of the dc model of a transistor network satisfies a certain dominance 
condition, then the network cannot be bistable. Therefore, a fundamental 
restriction on the G matrix of an interesting class of switching circuits is 
that it must violate the dominance condition. 

I. INTRODUCTION 

For each positive integer n let 5" denote that collection of mappings 
of the real n-dimensional Euclidean space E n onto itself, defined by: 
F t 5 n if and only if there exist, for i = 1, • • • , n, strictly monotone 
increasing functions /, mapping E 1 onto E 1 such that, for each x = 
(x, , • • ■ ,x n Y t E n , F(x) = (/,(*,), ■ • • , /„(*„))'• 

The main purpose of this paper is to report on some results concerning 

l 
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properties of the equation 

F(x) + Ax = B, (1) 

where A is an n X n matrix of real numbers, F maps E n into E", and 
B e E". In particular, a condition to be satisfied by A is given which is 
both necessary and sufficient to guarantee that for each FeSF" and each 
B e E n there exists a unique solution of equation (1). 

We also study the problem of obtaining bounds on the solution of 
equation (1). These bounds show that (if F t CF" and our condition on A 
is satisfied) the solution depends continuously on B. The bounds are 
often of use in computing the solution by standard iteration methods 
such as the Newton-Raphson method. By appealing to a theorem of 
R. S. Palais it is shown that the bounds can also be used to obtain a 
theorem essentially the same as, but somewhat weaker than, our 
principal result. 

Several results can be found in the literature which specify sufficient 
conditions for the existence of a unique solution of equation (1). For 
example, if A is positive semidefinite then a special case of a theorem 
of Ref. 1 guarantees the existence of a unique solution of equation (1) 
for all those F e ST which have the property that the slope of each /,- 
is bounded from above and below by positive constants, and for all 
B e E n . This theorem also specifies that a certain iteration scheme will 
always converge to the solution. 

A theorem of G. J. Minty 2 , when applied to equation (1), also implies 
essentially the same result. The boundedness condition on the slopes 
of the functions /,• is not required by Minty's theorem. On the other 
hand, Minty's theorem does not provide a procedure for computing the 
solution of equation (1). 

In Ref. 3 it is proved that a sufficient condition for the existence and 
uniqueness of a solution of equation (1) for all FeSF" and B t W is 
that A satisfy a weak row-sum dominance condition: 

a.-< ^ 2 I a a l> i = 1, "• , n* 

Other information concerning the location and the computation of 
the solution is also given in Ref. 3. 

The class of matrices satisfying the condition of our theorem (which 
is denned in Section III and denoted by P ) includes all positive semi- 
definite matrices as well as all matrices which satisfy any one of several 

* Appendix A contains a simpler proof of a similar result and a proof of a new 
related result. These results specify convergent algorithms for obtaining the 
solution. 



PROPERTIES OF DC EQUATIONS 6 

dominance conditions. Many other matrices are included in P \ and 
since the condition of our theorem is both a necessary and a sufficient 
one, we are assured that P is the largest class of matrices A for which 
equation (1) has a unique solution for all F e 5" and all B e E n . 

II. NONLINEAR NETWORKS 

Equation (1) is often encountered in the study of nonlinear electrical 
networks. In the case of networks containing only resistors (that is, 
linear resistors with nonnegative resistance), dependent and inde- 
pendent sources, and two-terminal nonlinear resistors that are described 
by functions in ff l (diodes, for example), this is rather obvious. 3 Even 
for networks which contain more general nonlinear devices, however, 
equation (1) can often provide a convenient characterization. For 
example, D. A. Calahan shows in his recent book that the transistor 
network of Fig. 1 may be described by the equation 

Cr /JtVr/tT -i\~ 

I..ie - 1) 

I c .{e" V,/kT ~ 1). 
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if the Ebers-Moll model is used to represent the transistor. (See pp. 
13ff of Ref. 4.) In this equation /„ , I e , , q, k, T, and V ee all represent 
fixed real parameters. It is quite trivial to apply the theory of this paper 
(in particular, Corollary 3 of Section IV) to Calahan's example and 
prove that this equation has a solution, the solution is unique, and 
the solution depends continuously on V e . . We also show how bounds 
on the solution can be obtained. 
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Fig. 1 — Biased transistor-stage. 
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More generally, it is frequently the case that networks which contain 
transistors, as well as the previously mentioned linear and nonlinear 
elements, may be described by the equation 

TF(x) + Ax = B. (2) 

In this case, a: is a vector whose components are the voltages across 
the nonlinear resistors and the transistor base-emitter and base-collector 
voltages. The n X n matrix A is the ?/-parameter matrix of the linear 
n-port network which is obtained by removing all nonlinear resistors 
and transistors and setting the value of each independent source to 
zero. The function TF(x) describes the behavior of the nonlinear re- 
sistors and the transistors. It happens that the matrix T is nonsingular; 
therefore equation (2) can be put into the form of equation (1). 

Networks which contains inductors and capacitors as well as the 
memoryless elements already mentioned are of course described by dif- 
ferential equations. Even the study of such networks, however, can 
often lead to the consideration of equations of the same type as equa- 
tion (1) . One usually finds the solution of such an equation is necessary, 
for example, when computing the solution of the differential equations 
by using some implicit numerical integration formula. 

The problem of determining the equilibrium states of the above- 
mentioned dynamic networks is one in which the consideration of 
equations of type (1) often arises in perhaps a more direct manner. 
In this regard, if it happens that equation (1) has a unique solution, 
then the network cannot possibly be bistable. 

When the determination of equilibrium states of a transistor net- 
work leads first to the consideration of equation (2) , then as a rather 
direct application of our existence and uniqueness theorem it follows 
that if the matrix A satisfies a weak column-sum dominance condition, 



a,-,- ^ £ | a,-,- |, i = 1, 



,n, 



then T~ l A t P and hence the network has exactly one equilibrium 
state. This result and related results which are proved in Section IV 
have the following interesting corollary: One cannot synthesize a 
bistable network which consists of resistors, inductors, capacitors, 
diodes, independent voltage and current sources, and one (Ebers-Moll 
modeled) transistor — or even an arbitrary number of (Ebers-Moll 
modeled) transistors with a common base connection. 

The authors feel that in many respects the main contributions of this 
paper are in the techniques used to prove the results. For this reason, we 
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have not chosen to summarize all of the results at the outset and relegate 
proofs to later sections. But rather, the results and the proofs will appear 
in the order in which they will best illustrate the techniques developed. 

III. MATRICES OF CLASSES P AND P 

The following notation will be used throughout the remainder of the 
paper: The origin in E" will be denoted by 0. If D is a diagonal matrix 
then D > (D ^ 0) means that each element of D on the main diagonal 
is positive (nonnegative). 

In Ref. 5 and Ref. 6 M. Fiedler and V. Ptak define the classes of 
matrices denoted by P and P . They in fact prove that the following 
properties of a square matrix A are equivalent: 

(i) All principal minors of A are positive. 

(ii) For each vector x j* there exists an index A; such that x k y k > 
where y = Ax. 

(Hi) For each vector a ^ there exists a diagonal matrix D x > such 
that the scalar product (Ax, D z x) > 0. 

(iv) For each vector x ^ there exists a diagonal matrix H z ^0 such 
that (Ax, H z x) > 0. 

(v) Every real eigenvalue of A, as well as of each principal submatrix 
of A, is positive. 

The class of all matrices satisfying one of the above conditions is de- 
noted by P. Fiedler and Ptak prove that the following properties of a 
square matrix A are also equivalent : 

(i) All principal minors of A are nonnegative. 

(ii) For each vector x ^ there exists an index k such that x k 9* and 
XkVk ^ where y = Ax. 

(Hi) For each vector x ^ there exists a diagonal matrix D z ^ such 
that (x, D z x) > and (Ax, D x x) ^ 0. 

(iv) Every real eigenvalue of A , as well as of each principal submatrix 
of A, is nonnegative. 

The class of all matrices satisfying one of the above conditions is de- 
noted by P • 

The following theorems follow directly from the above definitions. 

Theorem 1. If A t P then for every diagonal matrix A ^ (A > 0), 
A + A e P (A + A e P). 

Proof: Let x 9± 0. Then, since A t P , there exists an index A; such 
that x k * and x k (Ax) k ^ 0. Thus, x k (Ax + Ax) k ^ (>0). □ 
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In particular, Theorem 1 implies that if A z P and A ^ (A > 0) 
then det (A + A) ^ (>0). 

Theorem 2. If A e P then A' 1 z P. 

Proof: Suppose A & P. Let x 9^ Q be given and let y = A~ l x. y 9^ 6 
since A -1 is nonsingular. Thus, there exists a diagonal matrix D > 
such that (Ay, Dy) > 0, which implies (x, DA~ x x) > 0, or (Dx, A~ l x)>0, 
or (.A -1 :*:, Dx) > 0. That is, for every x ?± 6 there exists D > such that 
(A~*x, Dx) > 0. Hence A -1 z P. □ 

Because of the similarity of the definitions of the classes of matrices 
P and P , one might conjecture that this proposition is also true: If 
A £ P , and det A ^0, then A' 1 z P . This conjecture is in fact true. 
Interestingly enough, however, its proof is not obtained as one might at 
first suspect, by simply modifying the proof of Theorem 2. Moreover, the 
proof of this conjecture does not even seem to follow directly from any 
of the above definitions of P . Rather, upon making the trivial observa- 
tion that for every diagonal matrix D > 0, det (A~ x + D) = det (A -1 ) 
•det (Z) -1 -+- A) -det (D), the conjecture is easily seen to follow from the 
fact that det (D + A) ^ for every diagonal D > if and only if 
A t P . This fact is a direct corollary to the proof of Theorem 3. 

IV. EXISTENCE AND UNIQUENESS THEOREM 



The following theorem is the principal result of this paper. 

Theorem 3. If A is an n X n matrix then there exists a unique solution 
of equation (1) for each FtiJ" and for each B t E n if and only if A e P . 

Proof: (if) Let A t P , F z 5 n , md B t E". The solution of equation 
(1) is then unique (if it exists) since if x and y are both solutions then, 
using the strict monotonicity property of F, there exists a diagonal 
matrix D > such that F(x) - F(y) = D(x - y). But [D + A](x -y) = 
6 and, by Theorem 1, D + A is nonsingular. This means that x = y. 

We prove the existence of a solution of equation (1) by induction. For 
k = 1, • • • , n, let 



F k (x) = 



/.w 



./*(**). 



A k = 



a kk _ 



B k = 



&, 



.b k . 



Clearly, A k z P , F k z 5*, and B k z E k . Also, it is clear that there exists a 
unique solution of F x (x) + A t x = B x for each F x z 5 1 and for each B x z E 1 , 
and that this solution is a continuous function of foj . 
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Assume that there exists a unique solution of F k (x) +A k x = B k for each 
F k t 3=*, B k e E k , and that this solution depends continuously on any 
scalar parameter 77 upon which B k depends continuously. Let the 
matrices A k _ k+l and A k+ltk be defined by 



A k , k+1 = 

At+i.k = [flt+1.1 "■■ o*+i.*l- 

Then, for every real number x k+1 , the equation 

F k (x) + A k x + A k , M x M = B k (3) 

has a (unique) solution which is a continuous function of x k+ 1 and of 77. 
Let the components of this solution be denoted by x { = m<(g*+i , 77), 
for i = 1, • • • ,k, and define the vector M k (x k+1 , 77) by M k = {m x , • • • , 

m k y. 

We now prove that the function 

<p(x k+l , 77) = A k+lik M k (x k+l , 77) + a k+uk+1 x k+1 — b k+1 (rj) 

is monotone increasing in x k+1 : Let x l k+1 , x 2 k+1 c E with x\ +1 < xl +1 . 
Then, if M 1 (M 2 ) denotes the solution of equation (3) when z fc+1 = 
x l+i ( x l+i)> we have 

F k {M 2 ) - F k {M l ) + A k {M 2 - M 1 ) + A k . k+1 (x 2 k+1 - xl +l ) = d. 

Because of the strict monotonicity of the function F k , however, there 
exists a k X k diagonal matrix A > such that 



Hence, 



Thus, 



F k (M 2 ) - F k (M l ) = MM 2 - M 1 ). 



M 2 - M % = -[A + A k ]- l A k , M (xl +1 - xl +l ). 



<p(xl+i) — ¥>(*Li) = («t+i,ft+i - A k+lik [A + A k ] 1 A k , k+1 \(x 2 k+1 - xl +l ). 



But then, from the easily verified relation 



det 



flft + l.Jfc+l A k + lik [A 1 A k \ A k<k+i — 



A' 



0-0 



+ A k+1 



det (A + A k ) 
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and since 



A: 



+ A* + , 



^ (Theorem 1) , 



det (A + A k ) > 0, det 

_0 • • • 0. 

and x 2 k+1 — x l k+1 > 0, it follows that <p(x 2 k+l ) ^ v(zj +1 ). 

Now since <p is monotone increasing and, obviously, continuous in 
x k+l , it follows that the left side of the equation 



/*+i(**+i) + <p(x k +i) = 



(4) 



is a strictly monotone increasing function mapping E 1 onto E 1 , and hence 
equation (4) has a unique solution. If x'° +1 denotes this solution then 

n»,(a£+i) 

m k (x° k+l ) 

is the (unique) solution of 

F k+i (x) + A k+1 x = B k+1 . 

We must now prove that this solution is a continuous function of any 
scalar parameter rj upon which B k+l depends continuously. It suffices to 
prove that x k+1 depends continuously on r\ (see equation (3)). This may 
be done as follows: 

Let x° k+1 be the solution of equation (4) corresponding to 77 = 77 . 
That is, let 

1k + i(x° k+l ) + <p(x° k+1 , v°) = 0, 

and let € > be given. Since j k+l is a strictly monotone increasing map- 
ping of E 1 onto E 1 , so is /J+, , and hence /J+j is continuous. Hence, there 
exists 5' > such that if | /*+i(z°+i) — /t+iO^+i) I < &' then | x° k+1 — 
Xk+i I < «• Since <p is a continuous function of 77, there exists 8 > such 
that I 77 - 77 I < 5 implies | <p(x° k+l , 7, ) - v(x° h+1 ,v) \ < 8'.If\ r? - n\< 
8, and 



fk+i&k+i) + <p(Zk+i 1 v) — 0, 



then, 



PROPERTIES OF DC EQUATIONS 9 

U*i(xUi) - U+i(.X k+1 ) + <p(a£ + , , ij) - v(x k+1 , r,) 

= — [<p(xl+i , 7j°) - tp(x° k+1 , rff]. 
But since both / fc+ , and <p are monotone increasing in x k+l , 

(x° k+l - S* + i)[/* + l(4 + l) - /*«(**+!)] ^ 0- 

and 

(*2 + i - *i + i)M&-i , >?) " <P(x k+l , V)\ ^ 0. 
Therefore, 

| (x° k+1 — Zfc + i)[/*+i(a&n) — /ft+i(aifc4-i)] I 

^ | (.r° +1 - a;*+i)[^(»*+i , V°) ~ <p(x° k+1 , rj)] \. 
Now, if x° k+l = x k+x then of course | x° k+l - x k+l | < e. Otherwise, 

I U*i(x° k+1 ) - /»«&«) | g | „(*!!♦, , v°) - dA+i , i?) I- 

But then, 

I u+i(x° k+1 ) - n +i (x M ) i < s', 

and hence | x° k+l - x k+l | < e. Thus, x k+1 is a continuous function of 77. 

((»% 1/) Suppose A 4 P ■ K det A < then for sufficiently small 
f > 0, det (fJ -{-A) < 0. For sufficiently large f, however, 

det (r/ + A) = r-det (/ + j A) > 0. 

Thus, since det (f/ + A) is a continuous function of f, there is some 
value of f > such that det (f/ + A) =0. For this value of f let 
F(x) = fix. Clearly, for this choice of F t ff", equation (1) cannot have a 
unique solution. 

If det A ^ 0, but A has a negative principal minor, we can still find a 
diagonal matrix A > such that det (A + A) = 0; however, in this case 
A will not, in general, be simply the identity matrix multiplied by a 
positive constant f . 

For some positive integer k < n let A have n k X k principal minor 
which is negative and let 

A (1> = diag [5, , • • • , SJ. 

Since the determinant of A + A is not altered if any two rows and then 
the corresponding pair of columns are interchanged we may, without 
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loss of generality, assume that the matrix A is partitioned as 

Ai A 2 



A = 



_A 3 A 4 j 



where A, is a k X /c matrix with det A x < 0. Let £ > be chosen so small 
that det (£/ + A x ) < 0, and let 8, = • • ■ = 8 k = £. Now, if 5 h+1 = 
• • • = 5„ = f > 0, then 



det (A (1) + A) = det 



£/ + Ay A 2 

. A 3 tl + 4 4 J 



- f ~* • det 



g/ + A, 



A, 



(1) 



Thus, for f > chosen to be sufficiently large, det (A ( 1} + A) < 0. 
(det (A (1) + A) -» r""*-det (£7 + A,) < as f -» °o.) Now, if for 
r] > 0, A (2> = 77/, then it is clear that for rj chosen sufficiently large, 
det [A (2) + A] = 7j"-det (7 + (1/t?)A) > 0. Thus, if 

A(i) = eA (1) + (1 - e)A (2 \ 

it is clear, since det [A(0) + A] > and det [A(l) -f A] < and since 
det [A(e) + A] is a continuous function on ^ e ^ 1, that there is a 
value of e > (0 < e < 1) such that det [A(«) + A] = 0. For this value 
of e, A(e) > is the required diagonal matrix. □ 

Notice that our proof shows that if F e 5" and A e P , then the solu- 
tion of equation (1) depends continuously on any scalar parameter 
upon which B depends continuously. The arguments of Section V show, 
under these assumptions on F and A, that the operator (F + A)' 1 is 
in fact a continuous map of E n into itself. 

In the proof of Theorem 3 we see that the uniqueness of the solution 
follows simply from the hypotheses that each /,- is strictly monotone 
increasing and that A c P . The additional hypotheses that each f t 
is continuous and maps E l onto E l are not necessary (continuity of 
each /,• is not explicitly hypothesized, but follows from the "mono to- 
nicity" and "onto" hypotheses). Hence, we have: 

Corollary 1. If, jor i = 1, • • • ,n, S { is a subset of E 1 , and if S = S x X 
• • ■ X S n , and if F(x) = (/ifo), • • • , /»(&,,)) ', where each /,- maps E 1 
into E 1 and is strictly monotone increasing on S t , then if A t P and B e E n , 
there exists at most one solution of equation (1) in S. 
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We now prove another interesting corollary of Theorem 3. We first 
define some additional notation. 

For each positive integer n let S n denote the collection of all subsets of 
ET defined by: S z S n if and only if S = S 1 X ■ • • X S n where, for i = 
1, • • • , n, S i C E 1 and S* has the same cardinality as E\ For each 
5CS" we define the collection $"(£) of functions mapping S onto E n by: 
F z ff n (S) if and only if there exist, for i = 1, • • ■ , n, strictly monotone 
increasing functions /,• mapping S' onto E l such that for each x z S", 
Fix) = (MsO, ••• ,/»(*»))'■ 

Corollary 2. If A is an n X n matrix and the collection ff"(S) is non- 
empiy //lew //iere extste a unique solution of the equation 

F 1 (x) + AF 2 (a;) = B (5) 

/or eoc/i F! e 5 n (^)» ^2 £ 5 n (S), and each B z E" if and only if A z P . 
Proof: Since F 2 z 5 n (S), F; 1 : E n -> S exists and F, o F; 1 z £F n . Thus, 
there exists a unique solution of equation (5) if and only if there exists 
a unique solution of 

F 1 (F; 1 (y)) + Ay = B. □ 

As special cases of Corollary 2 we have: there exists a unique solution 
of each of the equations 

F x (aj) + AF 2 (x) = B, 

and 

x + AF(.r) = 5, 

for each F 1 , F 2 ,F t 5 n and each 5 z E" if and only if A z P . 

In Theorem 3 (and Corollary 2) the hypothesis that each of the func- 
tions /, is an onto mapping is quite necessary in order to guarantee the 
existence of a solution for each A z P . In the following example all of 
the hypotheses of Theorem 3 except this one are satisfied: 

e 11 + Xi — x 2 = 1 

e xt — x x + x 2 = —2. 

It is of course impossible for these equations to have a solution since, by 
adding both sides, we find that the solution would have to satisfy 

e" + e" = -1, 

which is absurd. 
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Even though the functions /,• are not "onto," it is still possible to 
specify sufficient conditions for the existence of a unique solution of 
equation (5) [and equation (1)] by strengthening the hypothesis on 
the matrix A — namely, by requiring that A z P. This is the essence 
of Corollary 3. We first require additional notation. 

With S" defined as above, we define, for each S z S n , the collection of 
functions 5o(S) mapping S into E n by: F z $l(S) if and only if there 
exist, for i = 1, • • • , n, monotone increasing f mictions /,• mapping 
S' onto a connected set in E 1 such that, for each x z S, F(x) = (ji(xi), • ■ ■ , 
/»(x„))'. When S = E n we denote 5%(S) by 5 n . 

Corollary 8. If A is an n X n matrix then there exists a unique solution 
of equation (5) for each F t z ff;(S), F 2 t 5 n (S), or F x e $ n (S), F 2 t JFJ(iS), 
and for each B zE",if A z P. 

Proof: If F 2 e ff n (S), F; 1 : E n -> S exists and F t o F; 1 t 5$ . Thus, in 
this case, there exists a unique solution of equation (5) if there exists a 
unique solution of 

F 1 (F; i (y)) + Ay = B. (6) 

Now, since A t P, it follows from the fact that the determinant of a 
matrix is a continuous function of each of its elements, that there is a 
matrix A* z P C Po and an e > 0, such that A = el + A*. Hence, 
equation (6) is equivalent to 

F(y) + A*y = B, (7) 

where we have defined 

F(y) - WW) + Jy. 

But, since F t o Fj 1 t 5FJ} and el t ff", it follows that F e $ n . Therefore, 
since A* z P , equation (7) and hence equation (6) and hence equation 
(5) have unique solutions. 

The case when F t z $ n (S) and F 2 t $o(S) can be reduced to the case 
just considered by making the simple observations that, in this case, 
equation (5) has a unique solution if 

A-'F^x) + F 2 (x) = A~ l B 

has a unique solution, and A z P implies A' 1 z P (Theorem 2). □ 

In Corollary 3 a sufficient condition is given for the existence of a 
unique solution to say equation (1) when the functions f t which specify 
F are not necessarily mappings onto E 1 . That the condition (A z P) is not 
necessary is easily demonstrated by the counterexample: Let F z CF§ and 
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B c E a ; then the equations 

f t (xi) — x* - fcj , and f 2 (x 2 ) + x x = b 2 
have a unique solution in spite of the fact that the matrix 

-l" 



A = 







4P- 



This is true because the function f 2 (ji(x 1 ) — &,) is obviously a continuous 
monotone increasing function of x t , and hence the left side of the equa- 
tion 

/ a (/i(*,) - 5.) + *i - h (8) 

is a strictly monotone increasing mapping of E 1 onto E 1 . Thus equa- 
tion (8) has a unique solution. 

V. BOUNDED SOLUTIONS AND RELATED PROBLEMS 

For many systems whose behavior is described by an equation having 
the form of equation (1), the vector B may be regarded as the sys- 
tem's input and the vector x may be regarded as the system's response, 
or output. Thus, if a sequence B l , B 2 , B 3 , • • ■ of input vectors for the 
system is given, the corresponding sequence x 1 , x 2 , x 3 , • • • of output 
vectors is specified by equation (1). An important property that such 
systems might have is that of producing a bounded sequence of output 
vectors for each bounded sequence of input vectors; that is, the property 
that whenever an input sequence B , B 2 , B 3 , • • • is contained in some 
bounded region of E", then the corresponding output sequence x 1 , x 2 , 
x 3 , • • • (exists and) also is contained in some bounded region of E n . 
By considering matrices A which are not members of P Q , it is easy to 
demonstrate that all equations having the form of equation (1) do not 
have this property. For example, if f(x) = x + e x (/ e 5 1 ), then the se- 
quence of solutions of the equation f(x) + ( — l)aj = b is unbounded, 
even though the sequence b = 1, §, |, • • • of inputs is bounded. The fact 
that one must resort to matrices A which are not in P , and the fact that 
by choosing any A $ P , an example of the above kind can be constructed 
by an appropriate choice of F t £F", follows from our next theorem 

Theorem 4- If A is an n X n matrix then A t P if and only ij for 
each F e 5" and each unbounded sequence of points x l , x 2 , x 3 , • • • in E n , 
the corresponding sequence B , B 2 , B 3 , ■ ■ • (B = F{x ) + Ax , k = 1, 
2, 3, • • •) is unbounded. 
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Proof: (if) If A 4 P then, as shown in the "only if" part of the proof 
of Theorem 3, there exists a diagonal matrix D > such that D + A is 
singular. Hence, there exists some point p c FT, p 9^ 6, such that Dp + 
Ap = 0. Let Pi , the j'-th component of p, be nonzero. Let the diagonal 
elements of the matrix D be denoted by di , • • • , d n and let the mapping 
F e "5 n be defined by 

/,.(*,) = \ diXi ' fOT **'"' 

Io\.t,- + e 1 ', for i = j. 

If p,- < let e = 1, if p t > let e = — 1. Consider the unbounded 
sequence x 1 , x 2 , x 3 , • • • defined by re* = k-e-p, for k = 1, 2, 3, • • • . 
The members of the corresponding sequence B 1 , B 2 , B 3 , • • • are B k = 
(0, • • • , 0, e ktv ', 0, • • • , 0)', k = 1, 2, 3, • • • , where the j-th element 
of each B k is nonzero. Since for k = 1, 2, 3, • • • , k e p,- < 0, the sequence 
B\ B\ B 3 , ■■■ is bounded. 

(only if) Our proof of the "only if" part of Theorem 4 consists of 
proving Theorem 5 which is referred to later for another purpose. □ 

Theorem 5. Let F = (fi(-), • ■ ■ , f„(-))' t 5", A t P , and, for i = 
1, • • • , n, a { ^ /3< be given. There exist, for i = 1, • ■ • , n, real numbers 
7,- ^ 5, such that for any B m (&! ,•••,&„)' e W with a t ^ 0, ^ fi { for 
i = 1, • • • , n, if x satisfies equation (1) then 7, ^ x { ^ 5,- /or * = 
1, ■ • • , n. 

Proof of Theorem 5: We first prove a useful lemma. 

Lemma 1. Let f be a strictly monotone increasing mapping of E 1 onto 
itself. Let x, b, a, /3 be real numbers such that xf(x) ^ xb with a ^ b ^ /3. 
Then y ^ x ^ 8, where 7 = min {/ _1 (a), 0} and 8 = max \f~ l (p), 0}. 

Proof: Let a ^ ^ /3 and define 7 = min {/"'(a), 0} and 8 = 
max {/"'O^.O}. Let x satisfy z/(:c) ^ so. Then x(f(x) - 6) ^ 0. Clearly, 
7 ^ ^ 5 and hence if x = then 7 ^ re ^ 3. Ifx>0 then /(z) ^ 
6:^/3 which implies x ^ f l (0) ^ 8 and hence 7 ^ < x ^ 5. If 
z < 0, then f(x) ^ ^ a which implies x ^ / _1 («) ^ 7 and hence 
y ^ x < ^ 8. D 

(Proof of Theorem 5) Since A tP there exists /c, e {1, • •• , »} such 
that ZtXilrc)^ ^ and hence, 

z*,o tl = rct./i.CxiJ + a; tl (i4.-c) 4l ^ a:*,/*,^*,). 
Thus, by Lemma 1, there exist y£ } =7" } (/*, , «*,) and «£> = «£'(/*, , /3*J 
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such that yff ^ x ki ^ tff. Now if F„_, denotes the mapping of E"' 1 
onto E"~ l defined by 

F-i ■ (/i(0, ••• ,/*.-i(-), /*.«(•), ••• ,/.(■))', 

if A n _, denotes the (n - 1) X (n - 1) matrix obtained from A by 
deleting the fci-st row and column (note that A n .i t P ), if 

and if 

B n _, = (6, , • • • , b tl _, , b fcl+ i , • • ■ , b„) , 

then 

F n - X (x) + A n _,x = B n . l - a n - x x ki * 

Since A n _, t P„ , there is a k 2 e {1, • • • , fc,-l, fci+1, •■• , n) such 
that x k ,(A n - l x) k , ^ and hence, as before, 

Xk,Q>k, — dk,.k,X kl ) ^ X k j k ,(Xk.)- 

But, if a kt . kl g 0, then 

a k , - a k „ kl yff ^ bjt, - Ofc,.*^*, ^ /3*. - a k ,, ki S k \\ 
and if a*,.*, > 0, then 

a*. - fli.,1,**" ^ b *. - a*..*.**. ^ 0*. - fl *..*.T*l > . 

Therefore, by Lemma 1, there is a yff = -y*" (/*. » "*. — A*..*!?*") 
and ill' = C(/*, . ft. " a t ,. kl Sff) such that <y«> ^ **. ^ C if 
«*,.*, ^ 0, and similarly for a k , kl > 0. 

The above process may be repeated successively until the n pairs 
of real numbers yff, Sff, (i = 1, • • • , n) have been obtained. Thus, for 
any given B with a, £ &< £ j8< far t — 1, • • • , n, the components of 
the solution x of equation (1) will be bounded by these pairs of numbers, 
provided it is known at each step which coordinate fc, to choose. The 
appropriate coordinate choice, however, will in general depend on the 
particular solution x which is associated with the given B. For different 
input vectors B the appropriate choice will in general be different. There- 
fore, in order to obtain bounds on x which are valid for all B with a, ^ 
bi t £ pi , (i = 1, • • • , n) we must consider each of the n! permutations of 
the coordinates {1, • • • , n) and, for each one, generate the set of bounds 
Wfft 8ff' *' = 1, • • • , w} for v = 1, • • • , n!. We then define 7. = 

* In this equation x is understood to be (ii , • • • , Xk,-t , z*-,+i , ■ ■ ■ 1 *n)'- 
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min M": v = 1, ••• , n\) and 5, = max {5J 10 : v ■■ 1, ••• , n!} for 
i = 1, • • • , n. Then, for each B with a, ^ 6, ^ /3, for i = 1, • • • , n, 
we have that 7, ^ a;,- ^ 5, for i = 1, • • • , n, since at least one of the 
sets of bounds {7^', d 1 ^: i = 1, • • • , n\ must always apply. □ 

If the matrix A of Theorem 5 satisfies a stronger condition than 
A t P (that is, if A satisfies a weak row-sum dominance condition), 



da 



= a2 I a a l» f° r * — 1] 



it is possible to use a method that requires much less computational 
effort than that of Theorem 5 to compute the vectors 7 and 5 whose 
components bound the corresponding components of the solution of 
equation (1). This method of computing the bounds, a straightforward 
generalization of an idea presented in Ref. 3, is explained in Appendix B. 

From Theorems 3 and 4 we now have the result: Every bounded 
input sequence B 1 , B 2 , B 3 , • • • is mapped by equation (1) into a bounded 
output sequence x l , x 2 , x 3 , • • • , for each F e £F", if and only if A t P . 

In the proof of Theorem 5, the number of real numbers 7*'', 8^ 
which must be computed, in order to determine bounds for x, is In X 
(n!). At the expense of obtaining poorer bounds it is easy to reduce 
this number to 2n 2 . Suppose we compute, at the first step, the 2n 
numbers 7, (1 \ *{", ■■■ , 7< n , C and set X x = min {y^, • • • , 7» 1> L 
Mi = max {81", ■ ■ ■ , 8™}. Then, for each B with a, ^ 6, ^ /3, for 
i = 1, • • • , n, one of the components of the corresponding x will be 
bounded by \ x (from below) and n r (from above). We next compute the 
In numbers y\ 2) = 7, (2, (/, , «, - p, a) ), S| 2) = «"'(/< , Pi ~ q\ l) ), where 
p,-° = max {a,,Xi , a if ni : j 7^ i\, g|- n = min {a,,X, , a,,Mi : 3 9* i\, and 
denote the smallest y i2) by X 2 and the largest 8 ,- 2> by n 2 . Then we have 
bounds which apply for two of the components of the x which 
corresponds to any B with a, ^ 6,- ^ /3, for i = 1, ■ • • , n. By computing 

7! 3) = 7| 3> (/< , «. " p|" - p|- 2) ), 5< 3) = «. (3) (/. , ft - «!" - <?, <2> ), etc., 
the above process may be continued to obtain the numbers X, , • • ■ , X„ , 
n, , • ■ • , n„ . Each component of the x corresponding to any B with 
<Xi S bi ^ /3, for i = 1, • ■ • , n will be bounded by X = min {X, , • • • , X n } 
(from below) and /* = max {jUi , ■ • • , n„] (from above). 

A matter that is closely related to the proofs of the above theorems 
on the boundedness of solutions of equation (1) is that of proving: For 
each F t 5" and each A t P the solution x of equation (1) is a continuous 
function of the vector B. It is obvious that it will suffice to prove that 
for each F t IF" with F(8) = 6, and for each A e P , the solution x 
of equation (1) is continuous in B at B = 6. We then note that if / 
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satisfies the hypotheses of Lemma 1 and, in addition, if /(0) = then, 
due to the continuity of f"\ for every e > there exists f > such that 
if or, /3 in Lemma 1 satisfy — f<a^&^/3<f then y, 8 in Lemma 1 
satisfy — e<7^.r^5<e. This observation may be used to incor- 
porate a simple "e-5 argument" into the steps of the previous paragraph 
to show that when F(6) = 6 then for arbitrary e > 0, one can determine 
f > such that \\B\\ < £ implies || x \\ < e. 

At this point we return to the matter of the existence and uniqueness 
of solutions of equation (1). We state first a theorem of R. S. Palais 
(Ref. 7 — see also the Appendix of Ref. 8) which shows the connection 
between the concepts of existence and uniqueness of solutions and the 
boundedness of solutions. 

Palais' Theorem. Let /i , • • • , /„ be n continuously differentiable real 
valued junctions of n real variables. Necessary and sufficient conditions 
that the mapping f : E* — » E* defined by f(x) = (fi(x), ■ • • , f„(x))' be 
a diffeomorphism oj E" onto itself are: 

(i) det [dfi/dXj] never vanishes, 
(ii) Iim|.i-No \\1(x) || = ». 

Palais' Theorem may be used to prove a result which is almost 
equivalent to our Theorem 3, that is: 

Theorem 6. If A is an n X n matrix then there exists a unique solution 
of equation (1) for each F = (/i(.r,), • • • , /„(x n ))' with continuously 
differentiable, strictly monotone increasing functions /,• which map E 
onto itself, and whose slopes are everywhere positive, and for each B c E", 
if and only if A t P . 

A proof of Theorem 6 which is independent of our Theorem 3 is 
easy to construct: For all A t P , the rather trivial Theorem 1 guarantees 
that condition (i) of Palais' Theorem is satisfied, and Theorem 5 
guarantees that condition (ii) is satisfied. If A 4 P then a choice of F 
such as is specified in the "if" part of the proof of Theorem 4 provides 
a case in which condition (ii) of Palais' Theorem is violated. 

VI. SUFFICIENT CONDITIONS FOR A t P OR P 

For a given matrix A, it is not in general an easy task to determine 
whether or not A satisfies any one of the four equivalent conditions 
of Fiedler and Ptak which are given in Section III and which serve to 
define the class of matrices P (or the conditions which define P). This 
is particularly true when the order of A is large. For this reason, we 
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now give several conditions which are sufficient to insure that a matrix 
A is in P or P (and which are not so difficult to verify). 

Suppose it were known that every eigenvalue of A as well as every 
eigenvalue of each principal submatrix of A had a nonnegative (positive) 
real part. Then this would guarantee that A t P Q (P). This is the main 
idea involved in the following theorem. 

Theorem 7. If any one of the following inequalities is satisfied by the 
elements a,, of the matrix A, for alii = 1, • • • , n, then A t P . 

(«) *« fe aS^E I *# D**. pfcl, P^ + ff" -1, 

a, positive numbers satisfying ^ (1 + «.-)"' = lj 

t-i 

(m) a,< ^ a max | a<, |, a positive satisfying 

E ( Z I <*.,- I (max | a u |) _1 } ^ a(l + a), (0/0 = 0). 

7/ any one of the above inequalities with ^ replaced by > is satisfied for 
i = 1, • • ■ , n, then A e P. 

Proof: If the right-hand side of any of the above inequalities is 
denoted by the nonnegative number r, then it is well known that all 
of the eigenvalues of the matrix A are contained in the union VJ { C,- : 
i = 1, • • • , n\ of the disks C,- = [z:\z- a,,- | ^ r,}. 9 But the condition 
a,-,- ^ (>) r,- guarantees that if z e C< then Re(«) ^ (>) 0. Thus, 
each of the eigenvalues of the matrix A has a nonnegative (positive) 
real part. The same is true of each eigenvalue of every principal sub- 
matrix of A, for if one of the above inequalities is satisfied by the 
elements of A it is also satisfied by the elements of any principal sub- 
matrix, n 

VII. COMPUTATION OF THE SOLUTION 

At present, the authors know of no single computational algorithm 
which is guaranteed to yield the solution of equation (1) for all F c JF", 
A e P , B t E". However, there are several ways that the solution may 
be computed for large classes of such equations. 

If, for example, the matrix A satisfies either a weak row-sum or weak 
column-sum dominance condition (inequality (i) of Theorem 7 with 
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either a = 1 or a = 0) and if F e $ B with, roughly speaking, the slopes 
of each /, bounded from below by some positive constant, then it can 
be shown (see Appendix A) that an algorithm for computing the solu- 
tion can be obtained by the use of Banach's contraction-mapping 
fixed point theorem. 

If the matrix A is positive semidefinite then, as mentioned in Sec- 
tion I, the existence of a unique solution of equation (1) for all Fcf 
follows from the earlier work of Sandberg and Minty. If, in addition, 
there exists for i = 1, • • • , n, positive constants a, and /3, such that 



«, ^ 



u — V 



for all u 9± v, then Sandberg's iteration scheme (also resulting from 
an application of the contraction-mapping fixed point theorem) can 
be used to compute the solution. 1 In this regard, if the techniques 
of Section V are first used to obtain bounds on the location of the 
solution then one could modify equation (1) by changing the nature 
of the functions /, outside the domain in which the solution is known 
to lie (but still keeping the /, strictly monotone increasing from E 1 
onto E l ) and obtain a new equation which has the same solution as 
the original equation. By doing this, the functions /,• in the new equa- 
tion might be made to satisfy the above inequalities in cases where 
this was impossible for the original /, . Also, even if these inequalities 
could be satisfied for the original equation, larger values of a { and 
smaller values of fi t might be used for the modified equation. This can 
result in a more rapidly converging iteration process (see Section VII 
of Ref. 3). Similarly, the bounds can be used to improve the perform- 
ance of other iteration schemes. 

In case A t P is not positive semidefinite, it might be that there exist 
diagonal matrices A, , A 2 > such that A,AA 2 is positive semidefinite. 
If such matrices can be found, then Sandberg's iteration scheme could 
be used to compute the solution of the equation 

A,F(A 2 :r) + A,4A a z = A,£, 

from which the solution of equation (1) may be obtained directly. 
Unfortunately, it is not the case that such A! , A 2 > exist for all 
A i P . For example, it is quite easily verified that for 

1 
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even though A t P , the matrix A,AA 2 is not positive semidefinite for 
any choice of Ax , A 2 > 0. 

It is easily verified, however, that appropriate A t , A 2 > can be 
found for all 2 X 2 matrices A e P except those for which 

(i) a n a 2 2 = 0, 

and 

(u) a l2 a 2 i = 0, 

and 

(Hi) either a, 2 ^ or a 2 , ^ 0. 

In particular, for all nonsingular 2X2 matrices A t P , appropriate 
A, , A 2 (A 2 = I) can be found. Thus, Sandberg's iteration scheme could 
be used, for example, to compute the solution of the example problem 
of Section II which was taken from Calahan's book. 

VIII. APPLICATION TO EQUATIONS FOR TRANSISTOR NETWORKS 

In this section some of the above theory is applied to the equations 
which describe the behavior of electrical networks containing tran- 
sistors. By the word transistor we refer to the three-terminal device 
whose equivalent circuit is shown in Fig. 2.* Considering the tran- 
sistor as a nonlinear two-port network, the following equations which 
express the port currents in terms of the port voltages follow imme- 
diately from inspection of Fig. 2: 



:h 



1 — a i2 



1M 



_—«21 1 _ -Ufa). 

We assume, as is the case for the usual large-signal model of a physical 
transistor, that < a 12 < 1, < a 2l < 1, and that both of the func- 
tions /x and / 2 are continuous and strictly monotone increasing. The 
character of the functions /, and f 2 which describe the transistor's 
nonlinear conductances will depend on whether the transistor is des- 
ignated as NPN or PNP. We shall, however, have no occasion to dis- 
tinguish between these two cases in what is to follow. 
Suppose an electrical network is synthesized by connecting together, 

* In some respects this equivalent circuit is an ideal model of a transistor. 
Nevertheless, since this model is often used in the design and the computer 
analysis of transistor networks, consideration of it is important. The presence of 
series resistance at the base, emitter, and collector terminals of a transistor will 
be considered by the authors in another paper. 
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Fig. 2 — The equivalent circuit of a transistor. 

in an arbitrary manner, any (finite) number of transistors, resistors 
(that is, linear resistors with nonnegative resistance) , voltage sources, 
current sources, and nonlinear resistors which are described by strictly 
monotone increasing conductance functions (and which we shall hence- 
forth refer to as "diodes") . Suppose the network contains n transistors 
and d diodes. For k = 1, ■ • •, n, let x 2k -i, Xzk, 2/2»-i> and y 2k denote 
the voltage and current variables v lt v 2 , ii, and i> , respectively, for 
the fc-th transistor. For k = 1, • • • , d, let x 2n+k and y 2 „+k denote the 
voltage across, and the current through, the fcth diode. Let these 
variables be related by y 2n + k = f2n + k (x 2n + k ). Then, if x = {x x , •••, 
x 2n+d ) * and y= (y lt • • -, y 2 „+a) ', we have 

y = TF(x), (9) 

where T = diag(7 T 1 , To), with T x a block diagonal matrix with ?! 2 X 2 
diagonal blocks of the form 

1 — «12 



■a 2 i 



1 



and T 2 a d X d identity matrix. The nonlinear function F has the 
form ^ 0c) = (fi(xi), ■•-,f2n+d{x 2n+ d)) t - 

Consider now the (2n + c/)-port network of resistors and independ- 
ent sources which is formed from the original network by removing the 
transistors and diodes. If the ^/-parameter matrix G of this (2n+d)- 
port exists then we have the additional equation relating the vectors 
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x and y : 

y = -Gx + u, (10) 

where u is some vector of constants which is, in general, nonzero since 
sources are present in the {2?i+d) -port. 

Combining equations (9) and (10) we obtain 



TF(x) + Gx m u. 



(11) 



Now T is a nonsingular matrix and hence, if equation (11) is multiplied 
by T~\ we obtain an equation having the form of equation (1). If 
the matrix T~ l G e P then, by Corollary 1, there exists at most one 
set of transistor and diode voltages satisfying equation (11). Moreover, 
if each of the nonlinear functions describing the transistors and diodes 
in our network maps E 1 onto E 1 , or if T~*G e P, then Theorem 3, or 
Corollary 3, guarantees the existence of a unique solution of equation 
(11). 

We have been careful to distinguish between the case when our 
theory guarantees only the uniqueness of a solution and the case 
when it guarantees both the solution's existence and its uniqueness 
for the following reason: In the analysis of transistor networks the 
nonlinear functions which are used to describe diodes or to describe 
the nonlinear conductances in the equivalent circuit of a transistor 
are often taken to be of the form 

f(x) = Z (e Xl - 1), 

where I and X are constants. The range of such a function is not the 
entire real line. Presumably, therefore, one can construct transistor 
networks having the property that if functions of the above type are used 
in a transistor's equivalent circuit then the network admits no solution. 
We now give a simple example of such a network. We wish to emphasize, 
though, that even for these networks whose equations may sometimes 
have no solution, our theory still guarantees that if T~ l G t P and if a 
solution of equation (11) exists, then it is unique. 

Consider the network of Fig. 3. For this network, equation (11) 
becomes 
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J>2, 
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Suppose a 12 = 0.5, a 2 i = 0.9, and g = 5.5 mhos. Then, the above equa- 
tion is equivalent to 
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Hence, v v and v 2 must satisfy 

/i(»,) + 5/ 2 (o a ) = 10(7. + 7 6 ). 

If we now assume that the transistor's nonlinear conductances are 
described by the functions 

fM - -I^'" ~ 1), 
MM = -Lie-*"' - 1), 

where the parameters I. , I c , X. , and \„ are each positive, then for 
all v, , v 2 we have 

fi(«,) +5/ 2 (y 2 ) < 7. + 57 e . 

Hence, if the values of the independent current sources of Fig. 3 are 
chosen such that 

7 U + h £ tV7, + \I C , 

then the equation for this network has no solution. 

Let us now consider the problem of determining whether or not, for 
a given network, the matrices T and G in equation (11) satisfy the 
condition T' l G e P (or T~ l G e P). (The existence of many transistor 
bistable circuits assures us that this condition is not always satisfied.) 



<£> 



o 



v 8 SjL c =fa(v 2 ) «2«t e (T) 



v,^i e = fi(v,) 



\^c\^) 



TRANSISTOR 



Fig. 3 — A transistor network whose equations may have no solution. 
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There is a large class of networks for which this condition is satisfied, 
and for which a simple inspection of the G matrix suffices to identify 
a member of the class. 

Since the matrix T satisfies a strong column-sum dominance condi- 
tion, that is, since 

Ui > T, I hi I for i = 1, ••• ,2n + d, 

the following theorem guarantees that if the matrix G also satisfies a 
strong column-sum dominance condition, then T~ l G e P, and that if 
the matrix G satisfies a weak column-sum dominance condition, 

On ^ Z) I 9a I. 

then T -1 (? e P and, hence, the above conclusions concerning the ex- 
istence and the uniqueness of a solution follow. 

Theorem 8. If the square matrix A satisfies a strong column-sum dom- 
inance condition and if the square matrix B satisfies a weak (strong) column- 
sum dominance condition, then A~ X B e P (P). 

Proof: Suppose A~ l B $ P . Then, by the main result of the "only if" 
part of the proof of Theorem 3, there exists some diagonal matrix 
D > such that det (D + A~'B) = 0. But det (D + A~ l B) = det (4 _I ) • 
det {AD + B), and det (A' 1 ) ^ 0. Likewise, det (AD -f- B) ^ 
since AD + B satisfies a strong column-sum dominance condition. 
Hence, A~ X B e P . 

With B strongly column-sum dominant, let 8 > be such that 
B — 8A also possesses the strong dominance property. Suppose that 
A~ l B — Si 4 Po . Then, as above, there is a D > such that A~ X B — 
hi + D = A~ l [B — 8A -f- AD] is singular, which is a contradiction. 
Therefore A~ l B - 81 t P , and, by Theorem I, A~ X B t P. □ 

IX. COMMON-BASE TRANSISTOR NETWORKS 

We now consider a special class of the networks which are comprised 
of transistors, resistors, diodes, and independent sources. We consider 
the class of all such networks for which there is a single node (called 
ground) to which the base terminal of each transistor is connected. 
Let us first consider a subclass of this class of networks; that is, let us 
temporarily assume that no diodes are present. For all networks in 
this subclass it is easily verified that when the G matrix for equation 
(11) exists, then it satisfies the above weak column-sum dominance 
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condition and hence, by Theorem 8, T~'G t P . This fact is made 
evident if we consider the network of resistors which is described by 
G (that is, the linear multiport to which the transistors are connected, 
with all sources removed) and first simplify this network by using the 
star-mesh transformation to remove all internal nodes. Of course for 
many networks of this sublcass G is strongly column-sum dominant, 
in which case T~ X G t P. 

It is clear that the networks for which the G matrix fails to exist 
are exactly those networks in which either one or more of the collec- 
tor or emitter terminals are connected, through the resistor network, 
directly to ground (that is, through a branch having infinite con- 
ductance) , or else two (or more) of the transistors' collector or emitter 
terminals arc connected directly together (through a branch of the 
resistor network having infinite conductance). These direct connec- 
tions can exist in the resistor network either because of corresponding 
short-circuits in the original linear multiport, or because of corre- 
sponding connections involving branches which contain only ideal 
voltage sources. 

If one assumes that each transistor in the network has a nonzero 
series resistance associated with both its emitter and its collector 
terminals (this assumption certainly being consistent with physical 
reality) then one need not be concerned about the possibility of the 
nonexistence of the G matrix since the situations mentioned in the 
previous paragraph cannot occur. We now show, however, that one 
need not rely upon this assumption in order to prove the uniqueness 
of the solution of the equations which describe the networks that we 
are considering. 

We have observed that the matrix G will not exist if and only if 
the linear multiport has fewer independent port voltages than it has 
ports. In this case we modify the nonlinear multiport in such a man- 
ner that we can break some of the connections to the linear multiport 
so that it then possesses a G matrix and hence can be described by 
an equation having the form of equation (10). The modifications 
to the nonlinear multiport which are called for are obviously the ad- 
dition of voltage sources between certain nodes, the values of these 
sources being the same as those of the voltage sources connecting the 
corresponding nodes in the linear multiport. This simple concept is 
illustrated in Fig. 4. Here, the network of Fig. 4a, containing a linear 
6-port, has been replaced by the "equivalent" network of Fig. 4b con- 
taining a linear 3-port, Although the G matrix of the 6-port does not 
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Fig. 4 — Example of a grounded-base transistor network. 



exist, it does exist for the 3-port which can be described by 
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i* 

We have shown that the above artifice allows an equation having the 
form of equation (10) to always be written to describe the linear 
multiport contained in our network. We now show that an equation 
like equation (9) can be written to describe the nonlinear part of our 



1 


-1 


Vi 




r ** 

-1 


1 


3 


Vi 


+ 


4 





1_ 


.08. 




. o. 



PROPERTIES OF DC EQUATIONS 



27 



modified network. The equation which we obtain is of the form y = 
PTF{P'x + C) with P an m X 2n matrix (m < 2n) and C a 2n-vector. 

Consider the equation which describes the nonlinear part of a com- 
mon-base transistor network before any of the above-mentioned 
modifications (that is, the addition of voltage sources) are made. 
This equation has the form of equation (9) with T being a 2n X 2n 
block diagonal matrix (recall that n is the number of transistors 
present) . Let us consider the effect on this equation of the modifica- 
tion of the network by adding voltage sources, one at a time. There 
are two different ways of adding voltage sources that must be con- 
sidered. 

Suppose a voltage source of voltage E is connected between nodes 
;' and k (with plus reference at node ;) , and suppose the connections 
between node ;' and the linear multiport are then open-circuited. This 
situation is illustrated in Fig. 5. Using the notation indicated in this 
figure, we have 



i' = TF(p), 



i y = i' 



for v t± j, k, 



ii = 0, 

i k = i'i + i'k , 

V; = v k + E. 

Let us now define the vectors v* and i* to be the (2n-l) -vectors 
obtained from v and i, respectively, by deleting the v, and t, elements. 
Then, if F*(v*) is the 2n-vector obtained from F(v) by replacing the 
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Fig. 5 — Typical modification of the nonlinear multiport network. 
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argument y, by v k -f- E, we then have that 

i* = T*F*(y*), (12) 

where the (2n— 1) X 2n matrix T 7 * is obtained from T by adding the 
j-ih row to the fc-th row and then deleting the ;-th row. Observe that 
T*F*(v*) can be written as QTFiQ'v* + R) in which the ;-th element 
of the 2n-vector R is E, all other elements of R are zero, and Q is obtained 
from the identity matrix of order 2n by adding the j'-th row to the fc-th 
row and then deleting the ;-th row. 

In case a voltage source of voltage E is connected between node ; 
and ground (with the plus reference at node j) and all connections 
between node j and the linear multiport are open-circuited, then we 
can again form equation (12) from equation (9) by simply replacing 
v,- by E wherever it appears in the argument of F, to form F*(v*), 
and deleting the j-th row of the matrix T, to form T*. In this case 
T*F*(v*) can be written as QTF(Q'v* + R) in which R is as defined 
earlier, but in this case Q is obtained from the identity matrix of order 
In by simply deleting the ;-th row. 

The above processes can be applied repeatedly to account for the 
addition of an arbitrary number of voltage sources to the nonlinear 
multiport. The resulting equation which describes the multiport will 
have the form 

y = Q P ■ • • Q 2 QrTF(Q[Q l 2 ■ ■ • Q' p x + C) 

m fP(x) 

with C some constant 2n-vector and each of the Q { obtained from the 
identity matrix of the appropriate order in one of the two ways described 
above. 

Consider equation (9) in which T is a square matrix. Due to the strict 
monotonicity of each component function of F, the mapping TF{x) 
has the following property: If p, q are arbitrary 2n-vectors then there 
is a diagonal matrix D > such that 

TF(p) - TF(q) = TD(p - q), (13) 

and furthermore, the matrix TD is strongly column-sum dominant 
(since T is strongly column-sum dominant). We now wish to show that 
a similar fact is true in the more general case. 

With m the number of rows of Q v , let p and q denote arbitrary 
?n-vectors. Then since there is a diagonal D > such that 

FWM ■■■ Qlp + C)- TOIQ; • • • Qlq + O = DQ[Q' 2 ■ ■ ■ Q' p (p - q), 
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we have 

TF(p) - TF(q) = Q p ■ ■ ■ Q 2 Q l TDQ[Q' 2 ■ • • Qj(p - g). 

The fact that Q„ ■ ■ ■ Q 2 QxTDQ{Q[ •■• Q' p is strongly column-sum 
dominant follows from the very easily verified proposition that the 
product Q k MQl (k = 1, 2, • • • , p) possesses that property whenever 
M does. 

Therefore if x l and x 2 denote two solutions of the "generalized equa- 
tion (11)," then [TD + GKs 1 - x 2 ) = d in which 

T = Q p • • • Q 2 Q,T and D = DQ[Q' 2 -■■(&. 

But TD, and hence TZ) + G, is strongly column-sum dominant and 
hence, nonsingular. This implies that x = x . 

We have now shown that in any network constructed from resistors, 
independent sources, and transistors having a common-base connection, 
the transistors' base-emitter and base-collector voltages are unique. It 
is a trivial matter to show that the same result applies when diodes 
are also allowed to be present in the network. 

Suppose the result was not true for some network containing at least 
one diode. Then there would be two different sets of voltages and cur- 
rents which satisfy Kirchoff's laws. Thus for each diode in the network 
there would be two (not necessarily distinct) pairs of points (v^ , i^), 
(v™ 7 i™) at which the diode is biased, corresponding to each solution. 
Letting / denote the strictly monotone increasing function which 
characterizes the diode we have i d Y) = f(v d iy ) and i d 2) = j(v d 2) ). But 
then, suppose the diode is replaced by the series combination of a 
resistor r and a voltage source E whose values are chosen so that the 
line i d = (l/r)v d - E/r passes through the points (v d l) , i ( d l) ) and {v d 2) , 
i ( d 2} ). (Due to the strict monotonicity of /, this can certainly be done 
with some positive choice of r.) Performing the above type of sub- 
stitution for each diode in the network, we obtain a new network of the 
type already considered. This new network would possess two different 
sets of transistor base-emitter and base-collector voltages (the same 
as before). This contradicts our previous result, and hence the previous 
result must apply, even when diodes are present in the network. 

To determine the equilibrium solutions of the differential equations 
which describe a network containing inductors and capacitors as well 
as the elements mentioned above, one must determine the solutions of 
a dc equation for a network of the above class. Thus, in summary, 
what we have shown is: One cannot synthesize a bistable network 
which consists of -resistors, inductors, capacitors, diodes, independent 



30 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1969 

voltage and current sources, and an arbitrary number of (Fig. 2) 
transistors having a common base connection (or, in particular, only 
one transistor) . 
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APPENDIX A 

Algorithms for Computing Solutions of Equation (1) 

In this appendix two algorithms for computing the solution of 
equation (1) are presented. It is proved that one of the algorithms will 
always converge to the solution of equation (1) if the matrix A satisfies 
either a weak row-sum or column-sum dominance condition (inequality 
(i) of Theorem 7 with either a = 1 or a = 0) and if, roughly speaking, 
the slopes of each /, are bounded from below by some positive constant. 
In each case the proof of convergence relies upon Banach's contraction- 
mapping fixed point theorem, and therefore also represents an inde- 
pendent proof of the existence and uniqueness of a solution of equation 
(1) for the conditions stated above. 

The following notation will be used: For fixed F z £F n , B z E n , let 
f(x) be F(x) — B\ also, if A is a given n X n matrix with elements a,-,- , 
we define the diagonal matrix D by D = diag [o n > (ha » * • • > fl nnL an d 
let A = A - D. 

Theorem A. If the n X n matrix A satisfies 

and if F e JT", B e FT, and if there exists some e > such that for each 
a, (3 zE 1 , c\ a — P\ ^ | /.(a) — /,(/3) \ for i = 1, • ■ • , n, then equation 
(1) possesses a unique solution, and if x° is an arbitrary point in E*, the 
sequence x°, x 1 , x 2 , • • • defined by 

x k+1 - (/ + D)-\-A)x k 
converges to the solution. 

Proof: Equation (1) may be rewritten as 

f(x) + Dx + As = 6. 
Hence, if the operation T: BT ->E"is defined by T = (/ + D)~\- A), 
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then the solution of the equation x = Tx is identical to the solution 
of equation (1). We now prove that the sequence x°, x\ x 2 , ■ ■ • converges 
to this solution by proving that T is a contraction. 

Let x and y be arbitrary points in E" and let g = Tx, h = Ty. Then, 
Kg) + Dg = -Ax and j(h) + Dh = -Ay. Thus, for i - 1, • •• , n, 

fM - h + Oitit = -(Ax),-, 
and 

/,(/*,) - b t +0,,/*, = -(Ay),. 
Subtracting, we obtain 

1M ~ U(hi) + a ti iSi ~ h) = (At/),- - (Ax).- . 
Since /,- is strictly monotone increasing, we have 

I UiSi) ~ U(h t ) | + a,-.- | g, - h< | = | (Ax),- - (Ay)< \, 
and hence, since e + a,, > 0, 

I* - hi i ^ — —- 1 <*■>* - m< i- 

Now, 

| (Ax),- - (Ay),- | = | 2 <*«(«* ~ 2/>) I 

^ (2 I a./ I) -max | x,- - y t |. 
!>• » 

Thus, denning the metric p on 2?" by p(x, y) = max | x,- - t/,- |, we 
have, for i = 1, • ■ • , n, 

But, since ^ 2 1 a <* I < a " + 6 » tnere exists X » ° - K < l ' suc ^ 
that | g, - /i. | ^ K-p(x, y) for i = 1, • • • , n, and in particular, p(Tx, 
Ty) = max | o, — h t | ^ K-p(x, y). Hence T is a contraction. □ 

Theorem B. 1 j then X n matrix A satisfies 

a,, ^ 2 I «»•■ li for * = !i *•• » n » 
and i'/ F t $ n , B t E", and i] there exists some e > swc/i //ia£ /or eac/i 
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a, jS c E 1 , e | a — | ^ | /<(<*) — /,(/3) | for i = 1, • • • , n, then equation 
(1) possesses a unique solution, and ij z° is an arbitrary point in 77", the 
sequence z , z 1 , z~, • • • defined by 

2 t+1 = _a(/ + 7))-V 

converges to some point z* and the solution of equation (1) is given by 

x* = (/ + DY l z*. 

Proof: As in Theorem A, the solution of equation (1) is also the solu- 
tion of x » (/ + D)~\-A)x. For each x e £"*, let z = (/ + Z))a: and 
hence a; = (/ + D)~ 1 z. Thus, x* is the solution of equation (1) if x* = 
(/ + D)~ x z*, where z* is the solution of z = — A(/ + Z)) _1 2. The theorem 
is thus proved if it is proved that the operator T = — A(/ + Z)) _1 is a 
contraction. 

Let P denote the operator (/ + D)~ , and let z and y be arbitrary 
points in E n . Then, proceeding as in the proof of Theorem A, we obtain 

| (Pxh - (Py) i | ^ — ~- I */ " Vi li for i = 1, • ■ • , n. 

Thus, if <7 = Tz and /i = T?/, then for i — 1, • • • , n, 

Qi = ~ X a„(Paj), and fc« - - 2 a ii {Py) i . 

Hence 

I g t - h, : | = | £ cCCft), - (Py),) | 

^ Z(l«w \-\(Px) i -(Py),\) 

s JS( | * ,l "^rn" l **" ft| ) - 

Therefore, 

= ±(e^ j -Ji*,-» ( i. 

i-i N Mi a u i e/ 

But, there exists /£, ^ K < 1, such that, for j = 1, • • • , n, 



lV) - a,-, + e 
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and hence 

£ | Qi - hi | £ K 2 I ** - 3/» I- 

i-i i=i 

Denning the metric p on S° by 

p(s, 2/) = Z I x t : - y. |, 
i-i 

we therefore have 

P (Tx, Ty) = _£ I 0. - A* I -a # -p(*> »)> 

i-l 

and hence T is a contraction. D 

APPENDIX B 

Determination of Bounds on the Solution of Equation (1) 

In this appendix we present a method for determining bounds on 
the solution of equation (1) when F t $ B , A is weakly row-sum dominant, 
and (for given a m ( ttl , • • • , «„)', = (ft , • • • , ft,)' t K) B = (6, , 

• • • , &„)' satisfies a, ^ 6,- ^ 0< for i = 1, • • • , n. The solution bounds 
are, in general, easier to compute than those of Theorem 5. The method 
presented here is a generalization of an idea presented in Ref. 3. 

The computation of the solution bounds proceeds in two steps. First, 
one solves each of the equations 

F(x) = a (14a) 

and 

F(x) = 0. (14b) 

Denoting the solutions of equations (14a) and (14b) by n = (mj , 

• • • , n n )' and v = ( Vl , ■•• , v„Y, respectively, and defining 

A = max{| Hi \, •■• , | /t« |, | Vi \, ••• , | vj}, 
and 

B' -(Zl-«|,-,S K>- 1)'. 

one then solves each of the equations 

F(x) + diag [a„ , • • • , a nn ]x = a - \B', (15a) 
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F(x) + diag [a„ , • • • , a nn ]x = + \B' . (15b) 

Denoting the solutions of equations (15a) and (15b) by y = (y { , 
• " f7n)' and 8 = (Si , ■ • • , 8 tt )\ respectively, one has y,- ^ .r° ^ 5< for 
i = 1, • • ■ , n, where .c° is the solution of equation (1) that corresponds 
to any B satisfying a { % 6, ^ & for i = 1, ■ • • , n. 

To prove that the components of the vectors y and 5, determined by 
the above procedure, are indeed bounds for the corresponding compo- 
nents of the solution x involves no more than a word-for-word repetition 
of the proof of Theorem 2 of Ref . 3, with several quite obvious modifica- 
tions. We omit the details. 
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